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ABSTRACT 


The equations of motion for turbulent two-dimensional 
magnetohydrodynamic flows are solved in the presence of finite 
viscosity and resistivity, for the case in which external forces 
(mechanical and/or magnetic) act on the fluid. The goal is to 
verify the existence of a magnetohydrodynamic dynamo effect which 
is represented mathematically by a substantial back-transfer of 
mean square vector potential to the longest allowed Fourier 
wavelengths. External forces consisting of a random part plus a 
fraction of the value at the previous time step are employed, after 
the manner of Lilly for the Navier-Stokes case. The regime explored 
is that for which the mechanical and magnetic Reynolds numbers are 
in the region of 100 to 1000. Ihe conclusions are that mechanical 
forcing terms alone c&nnct lead to dynamo action, hut that dynamo 
action can result from either magnetic forcing terms or from both 
mechanical and magnetic forcing terms simultaneously. Most real 
physical cases seem most accurately modelled by the third situation. 

®ie spatial resolution of the 32 x 32 calculation is not adequate 
to test accurately the predictions of the spectral power laws 
previously arrived at on the basis of the assumption of simultaneous 
cascades of energy and vector potential. Some speculations are offered 
concerning possible relations between turbulent cascades and the 
"disruptive instability." 


1 . 


JNTROrUTlON 


This is intended as the final paper in a series of three 
i -i n g with incompressible two-dimensional magnetohydroc’ynajni c 

turbulence. (A more exhaustive bibli ography has been given in the 
first two of these: Fyfe and Montgomery 1975; Fjr f^> Joyce, and 

Montgomery 19/6. ) In the first two papers, non- dissipative 
equilibrium predictions were algebraically derived and numerically 
tested, and an inverse vector potential cascade was predicted for 
the case of finite dissipation and external forcing. Our purpose 
here is to present the results of a numerical investigation of 
driven MHD turbulence for two-dimensional incompressible flow in the 
presence of finite conductivity and viscosity. 

Considerable interest in two-dimensional Navier-Stokes 
turbulence has been in evidence over the last several years. Much 
of the motivation has been meteorological, and stems from the 
prediction (Kraichnan 1967; Leith 1968; Batchelor 1969) of dual 
cascades : a simultaneous energy cascade to long wavelengths 

together with an enstrophy cascade to short wavelengths, both 
proceeding away from a source of localized excitation in wave number 
space. Kolmogoroff dimensional arguments lead to an omni-directional 
energy spectre. ~ in the enstrophy-caseading inertial subrange, j 

above which enstrophy is dissipated by viscosity at the same rate 

I 
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at which it is being supplied to the fluid. Ihe same Kolmogoroff 
arguments applied to the inverse energy-cascading inertial subrange 
lead to an omni-directional energy spectrum ~ k"'^, and an eventual 
accumulation of energy in the longest wavelengths allowed by the 
boundary conditions. Considerable meteorological support for this 
idea exists, and a number of numerical investigations of the 
viscous-decay initial-value problem have sought the k ^ subrange. 
However, to our knowledge, the only forced two-dimensional 
turbulence which has been simulated numerically, and for which it 
has been reasonable to expect an inverse cascade, is in the work of 
Lilly (1969, 1971). 

Lilly's approach to the problem has been followed rather 
closely in our generalizations to MHD. Lilly's work has recently 
been criticized by Herring et al^ (19?4) as spanning an insufficient 
range in wave number space to draw firm conclusions about values of 
spectral exponents in the two inertial subranges. We must concede 
the validity of this criticism. But while there is perhaps too 
much scatter in the computed spectra to draw firm conclusions about 
exponents, the existence of two different regions of the spectrum, 
the lowur one of which shows substantial back-transfer of energy in 
wave number space, seems to be beyond dispute in Lilly's 1969 paper. 
The two regions are also wholly compatible with the k"*^ and k”^ 
laws. In our opinion, the back transfer is the most physically 
significant effect in the entire problem. Since the analogous 


back-transfer of vector potential Vs what we most want to study in 
the MHD case, the inability to draw sharp conclusions about power 
law spectra seems a reasonable price to pay for its observation. 
Finally, Ully' s work is at this point still "the only game in 
town," and extending the range of allowed wave numbers to the point 
where the criticism of Herring et al. no longer applies would require 
a greater investment of computing capacity than anyone has as yet 
been willing to make. 

The essence of the criticism may be paraphrased as follows. 
Suppose the maximum allowed wave number is k and the minimum 

nicLx 

wave number is k . . Suppose the forcing terms added to the right 
min 

hand side of the Navier-Stokes equation are centered in the 
wave number space at a forcing wave number k^ Call the enstrophy 

3 1/6 

dissipation wave number k^ r k^ h (T)/v ) > where 7) = the average 

enstrophy per unit time supplied to the fluid, and v - the kinematic 

viscosity]. Then in order to have enough space between the various 

wave number ranges for both inertial subranges to be seen clearly, 

we require the following chain of inequalities: 

k . <c K. « k « h . (1) 

min F v max 

If each of these inequalities represents at least a factor of 10, 

x 

then the desired ratio of k /k is > 10 . For Lilly's Navier- 

max mm 

Stokes code (1969), this ratio was 32, and for our slower-running 

MHD code (a single run occupies about an hour of CDC 76OO time), 

k /k . is only 16. Numerical solutions for MHD turbulence in 
max' min 




which the basic grid has k nax A njin ~ 1000 lie beyond our present 

capacity. Navicr- Stokes computations with k /k . > 1000 are 

max' min 

presently at the threshold of numerical capability (Herring and 
0rs?.ag, private communication), but no computations have as yet 
been reported. 


It should be kept in mind throughout the following that in 

the MHD simulations under consideration, with two energy dissipation 

wave numbers (one constructed from viscosity, one from resistivity, 

in general unequal), no such chain of inequalities as (l) is 

possible either. Inevitably, the small size of our k will 

max 

affect the upper end of the calculated spectra in unphysical ways. 


and will prohibit our observing a proper dissipation range. 

In /a, we describe the general numerical procedure. In 
$ 3 - we ,HWress «» Nnvier-Stokes case (Lilly's problem); 

the bulk of the MHD results are presented in £ 4. Some spec ula tions 
on possible relations between MHD cascades and the "disruptive 
instability" are offered in ^ 5 . A general discussion of the 
results appears in ^ 6. 

It is well to stress the differences between these 


simulations and the earlier ones we have performed for the non- 
dissipative initial-value problem. Though •.hose are an essential 
preliminary to these, these are the realistic ones. For the non- 
dissipative equilibria, there is no energy transfer between 
different ranges of k at all, and the spectral levels never become 
independent of the finite k^, for which there is no physical 
determination. In the presence of forcing and dissipation, energy 
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transfer between various k values is the most central physical 
process, and though the system is continually trying to go towards 
an equilibrium state, it never succeeds. The interesting part of 
the problem is the balance struck between : (l) the input of 

excitations; (2) their dissipation; and (3) their thwarted attempt 
to reach equilibrium. 


'HP* 
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2. FORCED MAGNETOHYDRODYNAMICS 

As in our previous work, we employ periodie boundary conditions 
in the xy plane. Ihe MUD field variables have the following 
geometry : 

X = velocity field = (v , v , 0) 

x' y’ ' 

JB = magnetic field = (B , B , 0) = 7 x a 

x y »v 

a = vector potential = (0, 0, a ) 

z 

i = electric current density = (0, 0, j ) * 7 x B 

z ^ 

50 = vorticity density - ( 0 , 0 , o> ) = 7 y v . 

All quantities are expressed m terms of the appropriate dimension- 
less variables (Fyfe and Montgomery 1975). All quantities are 
independent of z, but are functions of x, y, and the time t. 
Supplementary constraints are that the v and B fields be solenoidal 

(7 y 0 and V • g _ 0 ) j these constraints are automatic when the 
equations of motion are written in the current-vorticity (J, „) 
representation. 

Assuming periodic boundary conditions in x and y, the Fourier- 
transformed MHD equations reduce to 
dao(k, t) . 

~ ^ E 6 ) 6 (£ + l » k) [®(r)»(g) - j(r)j( £ )] 
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and 


t) 

at 


= £) 


6 (£ + l “ W (j-Xfc) - «>(r)j(g)3 


- ^ 2 j(k, t) + g(k, t) 


(3) 


In Eqs. (2) and (3), all quantities have been Fourier-analyzed over 
a square box, so that, e.g. , 

t ) ex P (** ’ *) 

; z = ^ t ) ex P Ciis * *) •. 

®ie sums are over a large but finite set of k- values of the form 
k = 2tt (n , n )/L and L is the box dimension, n and n are integers, 
not both zero. 6 (g + r - k) is a Kronecker delta function of its 
argument, v is a dimensionless viscosity and n is a dimensionless 
resistivity. n/v is sometimes called the "magnetic Prandtl 
number." v” 1 , u” 1 are essentially Reynolds numbers, mechanical and 
magnetic, respectively. 

On the right hand sides of (2) and (3), f(k, t) and 
g(k, t) are the "forcing terms," mechanical and magnetic, 
respectively. They are regarded as the sources of the turbulent 
excitations which are dissipated by the viscosity and resistivity 
at the large values of k, after the excitations are shuttled through 


'pw wi iP‘-v i uwaiwpi jpvuu.p. 


“S™v 






T 




^-'v: 


'f **-> -r ■ 


"» - v --w \»««r *■*»> 


k- space by tlie nonlinear terms on the right hand sides of (2) and 
(5)- The coupling coefficients and M 2 are: 

**! (£» £) = • (rx E )(p‘ 2 - r" 2 ) 

^*p (r, p) = £ e • (r x p)|p + r| 2 p -2 r -2 . (1») 

The forcing terms f = f(k, t) and g = g(k, t) are regarded 
as given functions of k and c which model physical processes which 
may excite the MHD fluid. The range of physical processes which f 
and g may model is quite wide. For instance, they may be 
associated with a micro-instability of a non-MHD character which 
couples to the MHD motions as it grows. Also, f and g could 
represent stirring of an MHD fluid by a rigid conducting object or 
boundary. Or, f and g could be the result of an electric current 
driven in the fluid by an external transformer or capacitor bank. 

In the case of Navier-Stokes fluids, f can represent such different 
sources as a screen in a wind tunnel or the "baroclinic instability" 
(Pedlosky 1971) in meteorology. 

In any particular application, the functional form of f and 
g can be obviously quite complicated. However, in the spirit of 
homogeneous turbulence theory, we are seeking physical processes 
which are not specific to a unique situation but which are common 
to many different f*s and g'c. It is in this spirit that Lilly (1969) 
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gave a recipe for providing 
is relatively non- commit al, 


an f as a random function of time that 
and we follow his example for the MKD 


:ase. 


f(K, t) and g( 5 , t) arc set identically zero outside a fixed 

'ircular annulus ("forcing band") in V k y space. I* H & be the 
real or imaginary part of f or s for one of the allowed values of 
K in this band. The « k at time step n ♦ 1 is related to H, at time 

stet> n "by 


H. (n + 1) - 8 (n) +V 1 - B 2 J k (n) , 


( 5 ) 


rhcre r is a fixed real number between 0 and 1, and J, (n) is chosen 
)y a Gaussian random number generator with expectation value zero 
and variance such that the mean 0$ remains constant. In other 

words, the Fourier components are composed of a part which is a 
fraction of their value at the previous time step plus a random 
part chosen so that the expectation value of the magnitude of the 
forcing term is a constant. The evolution of the forcing field 
can be graphically portrayed as random walk over the surface of a 
hypersphere In the space of Fourier coefficients. The forcing term 
is programmed so that f and g can be chosen to be separately zero, 
to be uncorrelated, or to be identical. The limit corresponds 

to a temporally constant forcing field, and the limit R «= 0 to a 
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wholly random (and highly discontinuous) white noise field. The 
typical running value of R has been O. 95 . 

Ihe rest of the numerical scheme is very similar to that used 
by Fyfe, Joyce, and Montgomery (1976). It is a spectral-method code 
utiilzing the Fast Fourier Transform, along the lines of Orszag 
(1971) and Patterson and Orszag (1971)* The main accuracy checks 
have been tests of the conservation of the three non-dissipative 
invariants (Fyfe et al. 1976) with the dissipation and forcing terms 
turned off. To gain additional confidence in the program, we have 
chosen first to re-investigate the pure Navier-Stokes case 
[obtained by setting j s 0 and g h 0 in (2) and (5)] in erder to 
compare our results with those of Lilly ( 3969 ). This work is 
described in § 3* 

A final remark is in order concerning the way the numerical 
integration is carried out. Though the symmetry of Eqs. (2) and 
(3) renders them easiest to manipulate from a formal point of view, 
numerical calculations turn out to be easiest to perform on the 
equivalent pair of equations (Fouri<. r- transformed) 


S(W. 


Z = -v • 7a) + B • Vj + w 2 m + F 

~ z <v u 2 z 


at 


^ a z ? 

— , -V . 7a z < ^ a z + G , 


( 6 ) 


( 7 ) 
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where F, G are now the randan forcing functions. Because the 

equation for a , rather than j , is used, the Fourier transform of 
z z 

o 

G is not g(k, t), but rather g(k, t J/kS The implication is that 
if the forcing terms are to supply kinetic and magnetic energy at 
approximately the same rate, F and G are not to be comparable in 
magnitude, but rather the ratio of the Fourier transform of G to the 
Fourier transform of F must be of the order of the reciprocal of the 
square of the forcing wave number k^,. Unless specifically stated 
otherwise, it may be assumed in what follows that this has been 
taken into account and that when f and g are both non-zero, they are 
of the same order of magnitude. 

Ihe typical magnitude chosen for the forcing field is 

t 2 2 

£ (|f(k, t)j > = 16, where the summation is over those values of k 
between 55 and 70* A comparable value of £ (|g(k, t)J ) is chosen 
when magnetic forcing is present. Experiments have also been 
carried out with £ <|f(k, t)| ) = 64.0 and 4.0; the net effect seems 
to be to accelerate or decelerate the processes, and to raise or 
lower the limiting modal energies, but not to alter them 
qualitatively. 
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5. FORCED NAVIER- STOKES TURBULENCE IN TWO DIMENSIONS 

Three innovations not employed by Lilly (1969) have made it 

possible to improve some details of his computation, (l) We use the 

Qrszag- Patterson spt. tral method, not yet operative in 1969. 

(2) We average modal energies over allk-values corresponding to the 
2 

same k , instead of averaging around the perimeters of squares in 

k-space, as Lilly did. (3) We time average the modal energies once 

the total enstrophy has become approximately constant and once the 

rate of absorption of total energy has become approximately constant. 

All three improvements may reasonably be expected to reduce the 

errors in the spectrum slightly. 

Results for a typical Navier -Stokes run are shown in Figs. 1 

and 2. These refer to a forced, jurely Navier-Stokes run, starting 

with an empty spectrum. This run, as all others, is carried out on 

a 32 X 32 grid with k 2 =1, k 2 = 226, and a box size L = 2 tt. 

min max 

These parameters characterize all runs. The time step At is (128)”^, 

and this value characterizes all runs unless it is explicitly 

stated to the contrary. Some runs were carried out with 

At = (256) \ and there are some explicit comparisons described 

later for the MUD case. Fig. 1 shows the total computed enstrophy 

|«u(k, t)J 2 and computed energy 

« s k Ejr U( k > t )| 2 k " 2 = 4 Ej. t )| 2 versus time. 

it 
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The behavior of energy would be a linear increase with time and the 
behavior of enstrophy would be an approach to a constant, if the 
Kraichnan conjecture (1967) were given its most straightforward 
interpretation. 

The forcing wave number is approximately k^ = 8. More 
specifically, the forcing wave numbers are confined to a band 

p p 

extending from k = 55 to k = 70* Figs. 2a, 2b are time averaged 

modal energy plots. Figure 2a shows the time average of 

2 2 
Jv(k, t)| , averaged for given values of k , for 320 time steps 

2 2 

ending at t = 37* 5* Every value of k is plotted for 1 s k < 30, 

p 

every other point is plotted for 30 s k £ 100, and every fourth 

p 2 

point is plotted for k i 100. The allowed values of k get very 

dense at the higher values, and the above plotting convention is 

followed throughout the paper in order to keep the high-k modal 

energy plots from becoming unmanageably cluttered. Fig- 2b is a 

similar modal energy plot for an average of 1200 time steps ending 

with t - 75 . C. The straight lines drawn on the graph above and 

below the forcing wave numbers have slopes -4 and -8/3, respectively, 

2 -5 

corresponding to omni-directional energy spectra 2tt k < Jv(k) ( ) ~ k 

and k"'*/^. (Since the locations of the k , k points are far from 

x y 

isotropic at the lowest values of k, it seems misleading to convert 
the modal energy plots into plots of omni-directional spectra. ) 

I* 2 2 

Figure 2c is a plot of k |v(k)| vs. k , time averaged. 

Theoretically, this quantity should be flat in the ens trophy-cascading 


f 


1 (> 


subrange, and be proportional to k*/ 3 in the energy-cascading 
subrange. The straight lines have these slopes. 


Ibough it is apparent that the results are not inconsistent 
with the -3 and -5/3 predictions, the data do not permit the 
inference that these are the values. A least-squares fit of the 


expression Ah'" to k<|v(k)| 2 > for the data gives n . 1.49 for k 2 
below 55 and n - 2.07 for k 2 a 7 o. it has not been possible to observe 
the predicted build-up of the k 2 - 1 excitations above the predicted 
levelobtained from the k 3 ' 3 law. Though at the end of the run, 
the k 2 - 1 modes are still growing, the 9500 time steps represented 

in MLS ‘ 1 #nd 2 "* • lra *r <*«t 50* beyond the time over which 
the code will conserve energy and enstrophy below the 5* level for 

the inviscid unforced case. We believe that running longer would 
produce an arbitrarily large k 2 . 1 build-up, but the accuracy of 
the solution would have been lost. 


The enstrophy dissipation can be estimated as 
D - Svnk£ * kj, x [the energy dissipation] * 2.4, for the lMt 
1200 time steps. This gives a dissipation wave number 

* v * (2-VlO 9 ) 1/6 * 40, actually above k^. This makes the 

consistency of all but the highest-k points of Pig. 2 b with the 
-3 

’ * Pre<iit:ti0n qulte remarkable, and suggests that perhaps 
the inequalities (1) are too stringent. 

In looking at modal energy plots such as Pig. 2, on e should 
not be misled, by the logarithmic shrinkage of k space end by the 


plotting conventions, into underestimating the number of mooes in 

the upper reaches of the wave number space. For example, there are 

almost four times as many different k modes above k * 8 as there 

2 

axe below it. We have chosen tc plot k logarithmically, because 
in this particular application we are most interested in the back 
transfer to the lower values of k. It is the small-k part of the 
Fourier space that is emphasized by the log- log plots. 

2 

One should also keep in mind that each value of k has a 

2 3 

multiplicity of 2 to 8. That is, for each k , there are 2 to 0 
independent k values, and these degenerate k modes are always 
averaged over. 


I 
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4 . FORCED MHD TURBULENCE IN 1 WO DIMENSIONS 

It is well to remark at the outset the main qualitative 
difference that has been consistently observed between the MHD runs 
and the Navier-Stokes runs illustrated in $ 3- This is that the 
MHD case seems to be unavoidably "noisier ' than the Navier-Stokes 
case. As already noted (fyfe, Joyce, and Montgomery 1976) for the 
non-dissipative initial value problem, there is a slowing- down of 
magnetic energy transfer among the lower k's which makes time 
averaging less effective at smoothing out spectral fluctuations. 

Most of the MHD spectra are less smooth than those shown in Figs. 2 . 

A. The Case g = 0 (Mechanical Forcing ) 

The first case we consider is the case g(k, t) 2 0 . Biere is 
no significant amount of back transfer of vector potential for this 
case. The mean square vector potential can only decay. It is 
uncertain whether there is any physically realistic process which 
is satisfactorily modelled by setting g(k, t) 5 0 . For instance, 
mechanically stirring the fluid by constraining the velocity field 
to have a specified value at a solid surface also constrains the 
magnetic field at the same boundary (particularly if the solid has a 
finite conductivity) and so cannot be properly modelled without a 
magnetic forcing term. Nevertheless, we cannot entirely rule out 
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processes for which g 2 0, and we now present a few results for 
that case. 

Fig. 3 shows the enstrophy and energy vs. time for a run with 
g(k, t) s 0, v = = 10“ 3 , with a small initial "seed" magnetic 

field. Without magnetic field amplification, one might expect this 
case to yield something close to Navier-Stokes behavior, and it 
does. Energy increases linearly with time and enstrophy .fluctuates 
about a constant value. The corresponding modal energy plot (time 
averaged over 1200 time steps ending at t = 75*0) is shown in Fig. 4, 

this is to be compared with Fig. 2b. 

One interesting and nou yet quantitatively explained feature 

of this case is shown in Fi 0 . 5, a magnetic modal energy spectrum 
averaged over 1200 time steps ending at t = 75-0* (T* 1153 spectrum 

was flat, with all values = icf 6 at t = 0. Rrcm t = 0.0 to 
t = 75.0, the magnetic energy has decayed, by a factor of 2 X 10 . ) 
This steeply rising character for the very low level of initial 
magnetic excitations has appeared in all the purely mechanically 
forced runs we have done, including seme with still smaller seed 
fields. (The final ratio of magnetic to mechanical energy for this 
run is ~ 10~^, so that the errors in the computation of the 
velocity field are surely larger than the sizes of the magnetic 
components; but this is no reason to doubt the accuracy of the 
magnetic field calculation, assuming the velocity field to be given. ) 

A simple argument predicts qualitatively the growing 
|®(&)| 2 vs * behavior. Suppose wo regard the magnetic field as so 










*' * '* *m . -V • -s^ • "* -r 


S» -'■“'*■• :W* ^.^^P0^S«5Stf 


20 


* 


weak that it does not affect the velocity field v(x, t). Then the 
equation for the vector potential (omitting dissipation) 


aa z 

+ v • 7a = o 

ot ~ 2 
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shows that a^ is effectively convected like a passive scalar in a 

given turbulent velocity field, the effect of which will be to 

scramble the a field into progressively higher k modes. Batchelor 

(1959) has argued that under some circumstances, the omni-directional 

spectrum for such a passive scalar should approach ~ k - ^ behavior. 

This would imply an omni-directional magnetic field spectrum 
+1 

~ k . 3he actual rates of increase in such curves as Jig. 5 are 
typically higher than that, increasing more like k 2 ***. This 
discussion ignores the Alfven effect (discussed in Appendix A), 
which would predict a more rapid rise with k than that for a 
passively-convected a^. But we have no quantitative argument for 
the k * behavior shown in Fig. 5* The relation of the non - 
dissipative MHD and Navier-Stokes cases is discussed in Appendix B. 

B. The Case f - 0 (Magnetic Forcing \ 

The next case to be considered is the case of wholly 
magnetic forcing, f = 0. We are again uncertain as to the degree 
of realism to be assigned to such a model, but for the sake of 
completeness, we discuss it before passing on to the case where 
neither f nor g vanishes. Again, v = n = 10’^. 
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We begin from a completely empty spectrum 

<| g (*>| 2 - °> lx < S )| 2 - o ) • 

Fig. 6a shows the instantaneous magnetic nodal energies at t = 75*0. 
Fig. 6a reveals considerable scatter in the spectrum, but also 
significant back-transfer of magnetic energy to the longest wave 
length modes (recall that the forcing wave numbers are from 
k = 55 to k = 70). Fig. 6b shows the time-averaged (1200 time 
steps ending at t = 75*0) magnetic modal energies for the same run. 
The time averaging eliminates some of the fluctuations, but not as 

much as for the Navier-Stokes case. Die broken line has slope 

,-2/3 .X 

* > corresponding to an omni-directional energy spectrum ~ k 3 , 

though in no sense is it implied that the slope can be inferred from 

the points. Die corresponding kinetic modal energies are shown in 

Fig. 6c, time averaged over 1200 time steps ending at t = 75*0. 

Fig. 7 is a plot of total energy and vector potential vs. time for 

this run, and the cross -helicity P is shown also. As can easily 

be seen, vector potential is added at a constant rate, while the 

energy reaches a limiting value rather quickly and executes 

10$ fluctuations about that value. Diis is essentially what one 

expects for the dual-cascade behavior, and Fig. 7 should be compared 

with Fig. 3. Die cross-helicity is not identically zero, and 

fluctuates about zero with an amplitude of the order of .01. Dk: 




2 ? 

ratio of magnetic energy to kinetic energy at t * 75.0 is 1.74, so 
that the conversion of magnetic into mechanical energy has been 
substantial. 

By noting the rate of addition of squared vector potential 
(dA/dt ' 5 X 10' 4 from Fig. 7a) we may estimate the rate of addition 
of magnetic energy, which, since the energy added is only magnetic, 
enables us to estimate the rate of addition of energy, 

e(t)*kp dA/dt - .052. 

From this we may estimate the dissipation wave number 
k v = k - £e(t)/v 5 ] * (.32)* x io 2 * 75, 

above k 

max 

Another feature to be noted in Fig. 6a and Fig. 6b is the 
relative enhancement of the spectrum in the forcing band; the 
dissipation of energy away frem the forcing wave numbers is somewhat 
greater than it is for the Navier-Stokes case. 

A tentative explanation of this last-named effect is the 
following. The forcing wave numbers lie within a relatively narrow 
hand, and the first step in the excitation of the turbulence is the 
direct excitation of these wave numbers with k = ky by the forcing 
terms. The next step is the transfer of energy out of the forcing 

$ 

i 
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band by the nonlinear interactions measured by the coupling 
coefficients NL, ^ [Eq. (4)]. Now the ratio of the magnetic 
coupling coefficient to mechanical coupling coefficient for the 
interaction of two wave numbers r, g is 

e> Ie ♦ j | 2 

^ * p 2 . r 2 

For p, r about the sane size, this ratio is in general quite large. 
Urns for a narrouly forced band of excitations, magnetic excitations 
get out into the main body of the distribution faster than 
mechanical excitations. Since the excitations arriving at the 
higher values of k are dissipated first, the overall dissipation is 
expected to be higher for the case f = 0 than it is for the case 

g = 0. 


C. Mechanical and Magnetic Fomlnpr (r^Q i RtO) 

(toe surprise has occurred repeatedly, in the approximately 
12 runs (9600 time steps each) we have made for which we have used 
simultaneous magnetic and mechanical forcing: namely, a nearly flat 
spectrum of magnetic modal energies at all except the lowest values 
of k. This effect has persisted for different values of v, for 
correlated and uncorrelated f ang g; for variations in the 
magnitude of f and g; and for variations in the size of the time 
step. The effect is so pronounced as to he Identifiable as a 


?4 ... 

Qualitative difference Iran the Navicr-Stokes results (Figs. 2a, 2b, 
2c), from the mechanically' forced ME case (Figs. If, 5), and from the 
magnetically forced MID case (Mgs. 6a, b, c). Me are not confident 
■ that we have a correct explanation for the phenomenon, but we 
attempt one at the end of this section. 

Figs. 8 and 9 are not dissimilar to the cases previously 
discussed. Figs. 8 and 9 refer to the case of equal and uncorrelated 
f and g, for „■’. v = 10-->, and at = (lasf 1 . B,e forcing wave numbers 
•re again k 2 a 55 to 7 0. Fig. 8 is a plot of total energy, total 
cross helicity, and total mean square vector potential vs. time. It 
shows the energy e saturating and fluctuating about a saturated value, 
the cross helicity P executing small fluctuations about sero, and the 
mean square vector potential A increasing monotonically with time, 
much as in Figs. 7a, 7b. Fig. 9 shows the mechanical energy, the 
magnetic energy and their ratio R a [magnetic energy/kinetic energy] 

vs. time. Both energies saturate, and their ratio fluctuates about 
a value slightly larger than unity. 

mis behavior is so far consistent with the picture of the 
dual vector potential and total energy cascade we have been proposing. 
But a considerably more complex situation Presents itself when we 
turn to the spectral plots. Figs. 10a, 10b are time avenged metal 
energies vs. k 2 for the magnetic field and velocity field. 

Substantial dynamo action (back transfer of magnetic energy) ha. 
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obviously occurred, but except for the lowest magnetic modes, there 
is little structure of any kind, and comparisons with spectral 
indices cannot mean very much. 

Figs. 11a, b show exactly the same run, with the difference 

_2 

that now v - p = 10 . Increasing the dissipation by an order of 
magnitude is seen to have selectively dissipated the velocity field 
(there is no tendency for back transfer for Jv(k)j ), and to have 
substantially inhibited the back transfer of magnetic energy. The 
dynamo effect can, not surprisingly, be wiped out by making the 
viscosity and resistivity high enough. Neither e nor A show any 
systematic growth, for the situation depicted in Figs. 11, after 
about t = 5.0. 

In Figs. (11), the modal energies for the forcing wave 
numbers are off scale, at the level of > 1.4 x 10"^. This implies 
uhat by far the strongest nonlinear interaction is among the forcing 
band wave numbers, and these give the largest single contribution to 
the flux of energy into the high-k parts of the spectrum. We have 
previously noted that the ratio of the magnetic coupling coefficient 
to the mechanical coupling coefficient, 

M g (£» E> (E *lf 

«l<2> £> * p 2 . r 2 

is in general very large for similar-magnitude wave vectors £, r. 


*6 

!Biis implies that the highly excited forcing hand modes in Fig. 11 
are more effective at driving magnetic excitations outside the 
forcing band than mechanical ones, Ihe effect is strong enough to 
overcome the "Alfven effect" (Appendix A) and separate the magnetic 
and mechanical spectra as shown at the upper k values in Fig. 11. 

We have also done runs with uncorrelated f, g of equal 

o * 

magnitude and transport coefficients 4, v = 10 , 10 and 10 , 

-2 

10 . Perhaps unsurprisingly, the spectrum associated with the 

larger transport coefficient was depressed relative to the other, 
even at high k, by as much as an order of magnitude. By varying 
the 4/v ratio, the excitations can be made dominantly mechanical or 
dominantly magnetic. 

Examples of the kind of behavior exhibited in Figs. 8, 9, 

10 could be multiplied, but it seems redundant here to do so. We 
have a partial explanation to suggest for the behavior, though its 
confirmation will await simulations capable of considerably higher 
spatial resolution than k ^ = 16. We may estimate the "decay 

time" for wave number k by s (A)' 1 or (A) -1 , whichever is 
shorter, at the forcing wave number is about 15, ana at- the 

p 

maximum wave number (k mqy = 226), is about 5. We do not have an 
accurate expression for the rate of transfer of energy among 
different k values, but it may be faster than the eddy turnover 
time, r fc h (kjv(k)| J’ 1 , in any given region of k space. P*cm 
Fig. 10, is about 4 for the forcing wave numbers and about 2.5 


for *max* We have situation, then, that 5 for a?l k betwe 
kp and k^^. may reasonably be expected to measure th 'ime in 
which the non-dissipative equilibria are locally attained (Fyfe, 
Joyce, and Montgomery 19^6). The fact that the energy transfer 
among the allowed wave numbers may occur more rapidly than the 
m axi mum time of dissipation of any wave number could imply that for 
these simulations, the thermal equilibrium spectra (Fyfe and 
Montgomery 1975) are more relevant than the cascade spectra. What 
is being envisioned is a steady- state process in which energy is 
being added at a rate which is the same as the rate at which it is 
being dissipated, but that it rests in what is near to a quasi- 
stationary thermal equilibrium state. Biis picture is scmewhat 
substantiated by the dissipation spectra (|u>(k){ 2 , | j(k)| ? ) plotted 
in Fig. 12. These are seen not to have the characteristic behavior 
expected for cascade processes in their classical Kolmogoroff- 
Kraichnan form: a maximum, then decay near the upper end of the 

appreciably excited wave numbers. 

What the foregoing argument fails to explain are the 
differences in the MHD and Navier-Stokes dynamics. At the same rate 
of injection of energy, the Navier-Stokes fluid with identical 
transport coefficients is simply more efficient at disposing of the 
excitations at high k than the MHD fluid. The only explanation for 
this effect that ve have been able to arrive at is the already 
noted more rapid hlgb-k magnetic energy transfer. If correct 
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iae Implication of this for the accurate numerical simulation of 
full-blown dual cascades in MHD is somewhat depressing: namely, 
even a larger k^/k^ is required than for the Navier-Stokes case, 
Die current value of attainable ^ max A nln for the Navier-Stokes 
case of 1024 (Herring and Qrszag, private communication) is already 
on the threshold of present computational capacity. The MHD case 
may await the advent of the fifth-generation machines before it can 
be exhaustively studied. The more aware reader will have perceived 
uiat these considerations raise same questions about the adequacy 

of the spatial resolution of very nearly every MHD simulation that 
has as yet been performed. 

We close this section with some remarks on the effects of 
doubling the magnitudes of the forcing terms f, g. Pigs. 13 show 
modal energies for f, g twice as large, as in Pigs. 10. There is 
little change in the overall character of the plots, except for the 
enhancement of the k 2 = 1 magnetic excitation by a significant 

amount. Biis is the strongest "dynamo effect" we have observed in 
any single run. 

Other runs were carried out for the same parameters as 
Figs. 10 involve, but with a time step At = (256)’ 1 . No qualitative 
differences were observed in the spectra, though of course the 
spectra will not be identical since in general, different sets of 
random numbers will be involved in the forcing terms. 
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5. EXERCISE: A TURBULENT MHD CASCADE MODEL 

OF A "DISRUPHVE INSTABILITY” 


Under conditions of too high current or number density, 
tokamak discharges exhibit a behavior that has acquired the nsme 
"disruptive instability. " (For an encyclopaedic survey of the 
subject up to 1974, see Furth 1974; for recent experimental results, 
see, e.g., Hutchinson 1976 or Morton 1976; for recent theoretical 
approaches, see, e.g., Finn 1975 or Stix 1976). Essentially what 
happens is that the current channel in a cylindrical plasma column 
suddenly expands to fill the channel in a highly disordered way, 
in effect terminating the confinement (we ignore the toroidal aspects 
of the problem in this qualitative discussion). 

We wish to suggest here, without the extraordinary numerical 
effort that would be required in order to follow through on the 
details, a model, of this behavior which starts from a turbulent 
MHD cascade framework. Previously the phenomenon has been 
approached from a non- statistical MHD framework, usually (and not 
wholly successfully) from the point of view of linear MHD stability. 
The computational effort required to implement the program now to be 
out lin ed is not to be underestimated; numerical techniques would be 
required that cannot be described as "off the shelf" at present. 
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For computational convenience, virtually all the basic 
theoretical and computational studies of either Navier-Stokes or 
MHD turbulence have been carried out in rectangular geometry. 

Ihe turbulent fields have been expanded, assuming periodic boundary 
conditions, in the appropriate orthogonal functions, exp i k • x. 
Biough expanding in seme other set of eigenfunctions of the Laplacian 
would have increased the computational labor considerably, there is 
nothing fundamental about the restriction to rectangular geometry. 

For example, had we been interested in doing the present MHD 
calculation in the presence of rigid, perfectly-conducting circular 
boundaries instead of rectangular boundaries, an appropriate set of 
orthogonal functions in which to expand j and o> might have been 
exp(mKp), instead of exp ik • x [here (r, tp) are polar 
coordinates, m is any integer, and k^a is the nth zero or node 
of the Bessel function J ]. In fact, such a set of orthogonal 
functions would have had the advantage of permitting a non-zero 
value of the tolal current, JJ dx dy which in the case of 
periodic boundary conditions on the vector potential is constrained 
to be zero. (A non-zero value of this integral is a central 
component of tokamak confinement. ) Hie disadvantages would have 
been the loss of the utility of the fast Fourier transform and 
the Grszag-Ritterson spectral method: both severe losses. 

If the MHD equations, in the presence of other boundary 
conditions, are expanded in terms of another set of eigenfunctions, 
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their formal structure may be represented as 


dX. 

i 

dt 


jk 


C. X. 
ijk J 




- v. 
1 


X. 


+ F. 


( 8 ) 


Eq. (8) specialised to rectangular geometry and periodic boundary 
conditions, is just Eqs. (2) and (3). X^ is the ith element of a 
(possibly very long) column vector composed of all the independent 
expansion coefficients of all the necessary MHD variables. [In the 
problem considered in this paper, X^ represents symbolically all the 
o)(k) and j (k) for the allowed range of k. ] The are the coupling 

coefficients, given by Eqs. (4) for rectangular geometry; but in more 
complicated cases they would probably require numerical evaluation and 
storage. The is an external forcing ten-: (possibly random, 
possibly not). F^ might represent, e.g. , the electric field of the 
external transformer in the tokamak discharge; in that case, all but 
a few of the would vanish. Finally, v i is the viscous or 
resistive decay coefficient, the negative of the viscosity or 
resistivity times the corresponding eigenvalue of the Laplacian, in 
the present case. 

On the basis of the formal system (8), a scenario can be 
outlined which would lead to the qualitative behavior of the 
"disruptive instability." (A not totally dissimilar scenario can 
be imagined for the onset of turbulent flow of water down a pipe. ) 
Suppose that the electric field F^ is confined to one of the 
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low-lying values of i, say i^. By low-lying, we mean, corresponding 
to one of the smallest values of A solution of Eqs. (8) is then 

*io ” v io *io ( 9 ) 

with all the other X ± = 0 . In reality, the v i will have a temperature 

dependence, and are expected to be gradually reduced by the ongoing 

dissipation implied by Eq. (9). 3 he bigger F. is, the faster will 

10 

be the reduction. 

We have seen in the preceding pages that a turbulent cascade 
can be efficiently inhibited by increasing the dissipation rates. 

It may reasonably be inferred that a cascade may be launched, for a 
given forcing term ? iQ , when the fall below critical values. 
(Antecedents of this id°a appear in Kadomtsev 1975 and Sykes and 
Wesson 1976. ) 

For i ^ i Q , the initial stage of the developing turbulent 

cascade might be described by discarding all the terms in the sum 

in Eq. (8) except those involving X. . X. and v, could be 

10 10 i 

considered to a first approximation to be stationary in time. The 
growth of the other (in this phase and this phase only) would 
appear as a linear instability, for which the existence of a positive 
growth rate would involve a tlireshold condition involving , 

and the v^. As soon as the other X^ had grown out of the linear 
range, all the terms in Eq. (8) would have to be considered, and a 
full-blown turbulent cascade would be under way. Its accurate 
numerical description would be subject to all the stringent 
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requirements cn spatial resolution that our present simulation has 
been afflicted with, and without the shortcuts afforded by the fast 
Fourier transform and the Orszag spectral method- It is no 
exaggeration to say thrt no MUD simulation for any problem has yet 
been performed with spatial resolution small enough to avoid serious 
unphysical grid effects in a turbulent cascade* Whether this 
program, or its even more demanding three-dimensional analogue, could 
or should be carried out requires consideration of issues beyond the 
scope of this article. 

However, it is apparent that the above scenario does lead to 
a qualitatively correct behavior: excitation of the corresponding 

to coefficients of Bessel functions of higher index will manifest 
itself as a spreading of the j distribution over the whole cross- 
section. Hie enhanced transport that would surely result from the 
excitation of the vorticity coefficients would imply a transport 
of the KHD fluid to the walls at the same time. No obvious solution 
to the problem posed by the "disruptive instability" is suggested 
by the above considerations. Hie way to get rid of a turbulent 
MHD cascade is to make the plasma more resistive: cool it off or 

add more impurities. Neither message is particularly exhilarating 
in the context of controlled fusion research. 

Finally we remark that a somewhat more strained interpretation 
of the "disruptive instability" as an inverse cascade can be given. 

If the current distribution prior to the onset of the disruptive 
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Instability is not as smoothly centered as has been assumed, but 
rather is either closely confined to the center or is filamented 
(and the spatial resolution of the measurements of Hutchinson 1976 , 
for example-, do not seem to preclude this possibility, though they 
do not suggest it either), then the source of the excitations might 
be far out in Jt-spacc (or m, n space). ®en whe n a certain 
threshold was crossed, an inverse vector potential cascade to long 
wavelengths might result, or an inverse cascade of magnetic helicity 
for three-dimensional excitations (e.g., Fouquet et. al 1975 ). 
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6. Discussion 

It is considered that in the previous pages, it has been 
numerically established that dynamo action, in the form of 
substantial buck- transfer of vector potential, can occur for a two- 
dimensional magnetohydrodynamic flew with mechanical and magnetic 
Reynolds numbers large comparable to unity. Present computer 
capacity does not permit simulations on a large enough grid to be 
able to draw sharp conclusions about conjectured exponential 
power laws* The presence of external agencies, microscopic or 
macroscopic, has been modelled by the presence of forcing terms on 
the right hand sides of the equations of motion for the velocity 
field and magnetic field. It has been concluded that magnetic 
forcing terms must be present for the dynamo action to occur. Most 
processes which one might imagine as a source of turbulent 
excitations in an MKD fluid (such as friction with a solid surface 
of finite conductivity) seem best modelled by a forcing function 
involving both magnetic and mechanical components. 

Die translation of the above results into an experimentally 
testable format may be a slow and painful process. First of all, 
the diagnostic techniques for measuring fluctuating magnetic fields 
and velocity fields in even laboratory plasmas are primitive 
compared to those which have evolved for diagnosing Kavier-Stokes 


fluids (c.g. , Grant. Stewart, and Moilliet 19 &J; Fvcnkiel and 
Klebanoff 1967; Wiin-Nielsen 1967). For example, at present it 
appears to be difficult to measure the magnetic field interior 
to a tokamak discharge at all, let alone resolve it spatially over 
several orders of magnitude. Secondly, there are fonaidable 
psychological obstacles to the incorporation into plasma physics of 
the insights and techniques that have been achieved in the last 
fifteen years of fluid turbulence theory. "Plasma turbulence," as 
a subject, has evolved with a great deal of emphasis on linear 
instabilities and dispersion relations for normal modes of 
oscillation which are supposed to be a central component in the 
description of the turbulent field. None of these concepts plays any 
significant role in the modern theory of fully developed fluid 
turbulence. Readjusting the emphasis away from the preoccupation 
with linear apparatus (a twenty -five year enterprise that, in our 
opinion, has produced much of debatable scientific worth) 
may be slow indeed. Few industrial fluids, with the exception of 
those under storage in tanks or reservoirs, are stable in the sense 
which has been intended by plasma physicists; virtually all are 
turbulent. Fortunately, It is not necessary to have a good basic 
understanding of the turbulent flows involved in order to design 
fluid machinery. A proper and respectful distance has developed 
between the basic and engineering aspects of fluid mechanics* 

When a simile r division of plasma physics is perceived to be natural and 
desirable, the subject may be approaching maturity. 
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APPKNPIX A: VK2. ALFVEN EFFECT 


Hie Alfven effect war. first formulated by Kraichnan (l££) 
and ha., been invoked recently Pouquct ct au (1975) m edlJj ._ 


■tanped quasi-normal model computations for three dimensional MKD 
turbulent flows with helicily. n,e effect can be described more 
simply than has been done previously, end the description is 
relevant to seme of the other effects we have observed. 


Die effect concerns the- very high-k components of a 
turbulent MHD field. Ihe lower-k components are approximated by 
spatially-unifora velocity field and smgnetic field v 0 , & . Ihe 
Ihe highor-k v(k), B(k) are treated as small-amplitude perturbations 
on this spatially-uniform state. AU possible high-k fluctuations 
can then be shown to be a linear superposition of normal modes for 
Khich the time-averaged values of |v(k)| S and |§(k)| 2 are equal. 

In this approximation, one should see <|v(k)| 2 ) •_ <|B(k)| 2 > 
at high k, or equipartition of magnetic and kinetic energy. This 
has in fact been a prediction of the invlscid, perfectly-conducting 
equilibrium theory (fyfe and Montgomery 1975), end has been 
numerically verified for that case (We, Joyce, and Montgomery 
19 J 6 ). ihe equipartition at high wave numbers is a feature of 
Wga. 6, 10, and 15 of the present paper, and the fluid eppetrs to 
be striving for the cquiportition, even though it never achieves 
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it, in Pig. 5* However, we can see in Fig. 11 that the effect i: 
r.ot absolute, and the symmetry can be broken by increasing th*: 
dissipation, symmetrically or un cyrar. -‘trie ally , between the B and 
v fields. 

A second conclusion has been drawn about the Alfven effect: 

that the linear de-coupling of the fields v + B and v - B in the 

above approximation inhibits high-k energy transfer between wave 

numbers and thus flattens the energy cascade (an omni-directional 
_A J 

spectrum k 2 is supposed to replace the k behavior). About this 
conclusion ve are less sure. A transfer between the fields v + B 
and v - B may become more difficult, but a general uncorrelated 
non-helical forcing term supplies both fields at an approximately 
equal rate, and there seems to be no reason to suspect a decrease 
in the rate of transfer among different k’s for either field. Some 
of the present evidence favors an enhanced magnetic transfer. The 
conjecture must be left unresolved at present. It should be 
remarked that in such model calculations as the eddy-damped 
quasi-normal mold (Pouquet et al. , 1975 )> a range of exponents 
can be achieved by the choice of damping rate for the three- 
coefficient auto-correlations. The question of the implications 
of the Alfven effect for spectral laws must be decided on physical 
grounds outside the framework of models, since this choice is 
essentially arbitrary. 
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API‘£fn)IX B: gJUIIaBRIUM THEORY AT LOW A 

It was mentioned in Fyfe and Montgomery (1975) that the limit 
of the inviceid equilibrium theory as the magnetic excitation level 
approached zero (A -» 0) vis different from the inviscid Navier- 
Stokes equilibrium theory,, even though the dynamical equations 
approach each other in that limit. Here we digress from the mai n 
thrust of this paper to describe a numerical experiment which 
dramatically illustrates this subtlety involved in the limi ts. 

15ie experiment consists of finding the long-time modal energies 
for the velocity field for two inviseid, zero-resistivity initial- 
value problems. The veiocity-f.'eld coefficients are the same 
initially for the two runs. The only difference in the initial 
values for the two runs is the low level initial magnetic excitations 
in the second runj in the first run they are identically zero. 

For both runs, the initially non-vanishing velocity coefficients 
are for k = (3, 1), (3, 2), (2, 3), (3, -1), (3, -2), (2, -3), (1, 3), 
(1, -3). The initially non-vanishing vector potential coefficients 
for the second run correspond to k = (2, 2), (3, 0), (0, 3), (2, -2). 
Hie initial kinetic energy is h 5^ |v(k)| 2 = 2.94 for both runs, 
and the initial ragnetic field energy for the second run is 
® - O.O85, 2.8 9V» of the kinetic energy. The prediction 
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from the 1967 Kraichnan theory for the modal energy spectrum for 
the first velocity field is \|v(k)| ) = (a + pk ) with ot = -10.22 
and p = 11.17. 'jfce equilibrium theory for the second case predicts 
<|v(k)l 2 > = a" 1 , with a = 166.5* 

The numerically measured time averaged spectra at the ends 
of the runs are shown in Figs. l4a and 14 b. Bie theoretical 
curves are the solid lines. It is to us extraordinary that this 

addition of magnetic energy can, by destroying the constancy 
of the enstrophy, so alter the qualitative behavior of the 
velocity field. Indications are that any magnetic excitation, 
nowever small, can make this qualitative modification given 


sufficient time. 
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FIGURE CAPTIONS 

Fig. 1 Enstrophy fl(t) and energy e (t ) vs. time for pure Navier- 
Stokes case. A saturated constant ft(t) and a linearly 
increasing e (t) would characterize a perfect dual cascade, 
v = 10" 5 . 

2 2 

Pig. 2 Modal energies, | v(k) j vs. k , corresponding to the 

Navier-Stokes run shown in Fig. 1. Fig. 2a shows the 

average modal kinetic energies, averaged over 320 time 

steps ending with t « 37 . 5 . Fig. 2b shows modal kinetic 

energies averaged over 3.200 time steps (At = (128 ) -1 ) 

ending at t = 75*0. The forcing band is indicated by 
2 2 

arrows drawn at k =55 and k = 70 . Hie slopes of the 

broken lines correspond to omni-directional energy spectra 
-1-3 

~ k and k , respectively. Fig. 2c is a plot of the 

|i O 

"dissipation spectrum" k jv(k)| • The broken lines 

correspond to slopes for omni-directional energy spectra 
-I 1 -3 

~ k J and k , respectively. 

Fig. 3 Enstrophy fl(t) and energy e(t) vs. time for MHD run with 
mechanical forcing only (g = 0). Note the qualitative 
similarity to the pure Navier-Stokes case shown in Fig. 1. 

v ~ n = 10“ 5 . 


* 
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2 2 

Pig. 4 Modal energies |v(k}| vs. k for an MHD run with 

mechanical forcing only (seme situation as Fig. 3). The 
lines have slopes corresponding to omni-directional energy 
spectra - k~^ and k"^, respectively. The modal energies 
have been averaged over 1200 time steps ending at t * 75 * 0* 
Note the similarities to the pure Navier-Stokes run 
(Fig. 2). 

Fig. 5 Modal energies |B(k)| vs. k for an MHD run with 

mechanical forcing only (same run as Figs. 3 and 4 ) and a 
small initial seed field jB(k)| 2 - 10 all k. The 
modal energies have been time averaged over 1200 time steps 
ending at t = 75*0. This evolution to a sharply rising 
magnetic spectrum is characteristic of all runs of this 
type. 

Fig. 6(a) Modal energies |B(k)j 2 vs. k 2 f or a p, = v = 10 -5 MHD run 
with magnetic forcing only (f = 0). These values are 
instantaneous values at t ” 75*0# and their comparison 
with the time averages in Fig. 6(b) shows a typical level 
of fluctuation-reduction achieved by the time averaging. 
The broken line drawn below the forcing band has a slope 
corresponding to an omni-directional magnetic energy 
spectrum ~ k . 
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(b) Modal energies J B ( K ) | vs. k for the same MHD run shown 
in Fig. 6(a), but time averaged over 1200 time steps 
ending at t = 75.0 The broken line corresponds to an 

omni-directional magnetic energy spectrum ~ k"^. 

2 2 

(c ) Mortal energies |v(k)J vs. k corresponding to Fig. 6(b), 
again time averaged over 1200 time steps ending at 

t = 75.0. 

Fig* 7( a ) Total energy e(t) vs. time for the MHD run (f = 0) 

shown in Fig. 6. For an energy cascade, this number would 
saturate and become approximately constant in time. 

(b) Cross -helicity P versus time. 

(c) Mean square vector potential A(t) vs. time corresponding to 
the MHD run (f = 0) shown in Fig. 6. A pure inverse 
vector potential cascade would appear as a monotonic 
increase of A with time. 

Fig. 8 Total energy e(t), cross helicity P(t), and mean square 
vector potential A(t) vs. time for an MHD run with 
uncorrelated mechanical and magnetic forcing 
(f 0, g ^ 0) of approximately equal magnitude. For a 
dual-cascade situation in which no cross-helicity was 
being added, c(t ) would saturate, P(t) would fluctuate 
close to zero, and A(t) would increase linearly. 

H = v = 10 . 
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Fie* 9 Total mechanical energy, magnetic energy and ratio R vs. 

time for the forced MHD run shown in Figs. 8 and 10. 

Fig. 10 (o ) Modal energies | B(h)| 2 and (h) modal energies |v(k)| 

for the forced (f t 0, g t 0) MHD run shown in Figs. 8 and 
9. Substantial back-transfer of magnetic excitation into 
the loneest wavelength modes has occurred, but the 
fluctuations are larger than for either the f = 0 or 
g = 0 cases. Both (a) and (b) represent time averages 
over 720 time steps ending at t = 71.0. 

Fig. 11(a) Modal energies |B(k)| 2 and ft) modal ener « ies lxfe)| 2 vs * 
k 2 for a forced (f / 0, g i 0) MHD run which is the same 

as Figs. 8, 9, 10, except that the dissipation has been 
increased ton = v = 10" 2 Averages have been taken over 
1200 time steps ending at t = 62.25* The forcing band 
modal energies are off scale, being >1*4 * 1° • 
most unusual feature of Fig. 11 is the separation of the 
mechanical and magnetic energies at the highest k values. 
(Neither e nor A increases for this case after about 
t - 3*0.) 

Fig. 12 Dissipation spectra, (a) k 2 |v(k)l 2 an* ft) k ISOsM » 

for the MHD run (f /0, g i 0) shown in Figs. 8, 9, 10* 

The pure double-cascade theory would predict two straight 

lines for |B(k)| 2 vs. k 2 , the one below the forcing band 

with a slope of § and the one above the forcing wave number 

2 2 

with a slope of either 1/6 or 1/4, on a £n k |g(k)| ' vs. 


4f» 


to k' plot. Neither region is fit very well by these 
straight lines. If H m<|v /k were much larger, this 
spectrum should eventually reach a maximum and decrease 

p 

sharply with increasing k'. 

Fig. 13 (a ) Magnetic modal energies, a..l (b) kinetic nodal energies 
2 

vs. k for a forced MHli run (f / 0, g 4 C) similar tc 

Figs. 8, 9> 10, except that the magnitudes of f snti g 

have been doubled. Both (a) and (b ) ai time averages 

over 1200 time steps ending at t * 75.0. Fig. (13a) 

shows the greatest back-transfer of magnetic energy to 
2 

the k = 1 modes cf any of the r uns . 

Fig. lU (a) Modal energies, time averaged over 1200 time steps ending 
at t = 37*5 for pure Navier-Stokes initial- value problem 
with zero viscosity. The theoretical prediction is the 
solid line, ihe fractional variation over the duration of 
the run was 0. 4> for the energy and 5.4$ for the enstrophy. 
(b) Modal energies for the velocity field for the same 
situation as Fig. (14a), but with an approximately 3$ 
addition of magnetic excitation (u = v ® 0). The run 
is time averaged over 1024 time steps ending at t * 40.0. 
lhe theoretical curve is the solid line (flat, in this 
case). 
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